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pairs of Higgs doublets. We study the possibility of spontaneous CP violation 
in these scenarios and present a model where the origin of CP violation is 
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lO . We consider supersymmetric extensions of the standard model with two 
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lO ' soft, with all the complex phases in the Lagrangian derived from complex 

^ ■ masses and vacuum expectation values (VEVs) of the Higgs fields. The main 

p • ingredient of the model is an approximate global symmetry, which determines 

{iLii the order of magnitude of Yukawa couplings and scalar VEVs. We assume 

\ that the terms violating this symmetry are suppressed by powers of the small 

^ \ parameter epg = 0{mh/mt). The tree-level flavor changing interactions are 

^ , small due to a combination of this global symmetry and a flavor symmetry, 

but they can be the dominant source of CP violation. All CP- violating 
effects occur at order ej>Q as the result of exchange of a/most-decoupled extra 
Higgs bosons and/or through the usual mechanisms with an almost-ieal CKM 
matrix. On dimensional grounds, the model gives ex ~ Cpg and predicts for 
the neutron electric dipole moment (and possibly also for e'^) a suppression 
of order epg with respect to the values obtained in standard and minimal 
supersymmetric scenarios. The predicted CP asymmetries in B decays are 
generically too small to be seen in the near future. The mass of the lightest 
neutral scalar, the strong CP problem, and possible contributions to the Z 
decay into b quarks (the Rb puzzle) are also briefly addressed in the framework 
of this model. 



1 Introduction 



Although the standard model is today in impressive agreement with all particle physics 
data, its scalar sector has not been proven yet. A scalar sector defined by elementary 
fields seems to contradict the possibility of two very different mass scales (namely, the 
electroweak and the grand unification or Planck scales). Supersymmetry (SUSY) 
would offer an explanation for the stability at the quantum level of the different scales 
of the theory, provided that it is broken only by soft terms below the TeV region. A 
lot of attention has been paid to the minimal SUSY extension of the standard model 
(MSSM), which presents appealing features such as a consistent grand unification of the 
gauge couplings or a candidate for the cold dark matter of the universe. Experimentally, 
the MSSM has been so far flexible enough to avoid confiict with any measurement, but 
its most compelling prediction, the presence of a light neutral Higgs (lighter than the Z 
boson at the tree level), is still missing. 

The MSSM would also offer distinctive predictions for CP violating processes. Arbi- 
trary complex phases ip in soft gaugino masses and scalar trilinears would give fermion 
electric dipole moments (EDMs) well above their present experimental limits. This im- 
plies generically ip < 10^^ [@], a somewhat unnaturally small number. The MSSM 
would also predict unsupressed flavor changing neutral currents (FCNCs) unless there is 
some degree of degeneracy between squark masses (something which occurs for super- 
gravity Lagrangians with canonical kinetic terms) and correlation between the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix and its equivalent in the squark sector. In the usual 
SUSY scenario CP violation in K and B physics depends essentially on only one 
phase (the CKM phase 0), whereas the set of small phases in soft terms (uncorrelated 
to the family structure) may have experimental relevance only in fermion EDMs. This 
scenario, however, holds only for highly degenerated squark masses. In general, taking 
the experimental limit — M f^o^ K — K mixing one obtains that acceptable 
complex phases in gaugino masses may have an impact on the K system. For complex 
gluino masses this was shown in 0, and a model with small phases has been recently 
proposed in p. It was also shown ^ that due to large top-quark effects acceptable 
complex phases in chargino mass terms may also contribute to CP asymetries in the K 
system. 

A different approach to the origin of CP violation which is specially appealing in 
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SUSY models is the idea of spontaneous CP violation (SCPV) 0. All the phases in 
the Lagrangian (initially CP-conserving) would have their origin in a small number of 
complex vacuum expectation values (VEVs) of scalar fields. Moreover, the sizes of these 
phases could be correlated by approximate symmetries suppressing some couplings in 
the Lagrangian |^. Unfortunately, in the minimal Higgs sector of the MSSM there is no 
room for SCPV [rU, TT|, and hard CP violation is required. The possibility of SCPV has 
been also studied in SUSY extensions with singlet fields 0, where it can be obtained 
but seems to require certain ammount of fine tuning. 

The SUSY models with more than two Higgs doublets are an obvious extension of 
the MSSM ||13|. They are minimal in the sense that no new species are introduced, but 
just repeated. From the model building point of view there is no compelling reason to 
disregard them, and they could appear naturally in models with fermion-Higgs unification 
(like Eq) or left-right symmetric scenarios, where two bidoublets are required in order 
to obtain realistic fermion masses and mixings. Four Higgs doublet (4HD) models require 
an intermediate scale to be consistent with grand unification, but even this could be more 
in line with recent data on as{Mz) than the desert scenario |]T5|. Since more than one 
Higgs doublet couples to quarks of a given charge, a possible concern in this type of models 
is the presence of FCNC at the tree level. The experimental limits, however, can be easily 
avoided just by invoking the action of an approximate fiavor symmetry (see next section). 
On the other hand, a nonminimal scalar sector opens the possibility of SCPV and, in 
general, widens the parameter space relevant in low-energy precision measurements (this 
could be convenient, for example, if the anomalous value of Ri, persists). 

In this paper we present a 4HD SUSY model which seems to contain satisfactory 
answers to many phenomenological questions. CP violation appears softly, in complex 
Higgs masses and VEVs. The main ingredient of the model is a Peccei-Quinn like ap- 
proximate symmetry which determines the order of magnitude of Yukawa couplings and 
scalar VEVs. We define this symmetry in such a way that the additional pair of doublets 
has small VEVs with order one complex phases and is weakly coupled to all matter fields. 
As a consequence, the ratio mb/rrit, CP- violating effects in K physics, and the neutron 
EDM will appear suppressed by powers of the small parameter epg that parametrizes 
the violation of this symmetry. The CP asymmetries in B decays are predicted to be 
typically two or three orders of magnitude smaller than in CKM scenarios (the CKM 
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matrix in the model is essentially real), a signal that can be used to discriminate this 
4HD model with respect to the MSSM or the standard model. 

The plan of the paper is as follows. In Section 2 we write the generic Lagrangian for 
4HD SUSY models and review previous results on spontaneous CP violation. We show 
that a reahstic scenario for soft CP violation requires complex Higgs masses in the initial 
effective model. In Section 3 we define our model and minimize the Higgs potential. We 
show that the order of magnitude of Yukawa couplings, complex scalar VEVs, and the 
CKM complex phase are correlated by the approximate global symmetry. In Section 4 
we explore the implications of the model on K and B physics as well as on the neutron 
EDM. In Section 5 we discuss the resulting spectrum of scalar fields (in particular, the 
mass of the lightest neutral mode) and other possible phenomenological impacts of the 
4HD model. Section 6 is devoted to conclusions. Details about the minimization of the 
scalar potential can be found in the Appendix. 

2 Complex VEVs in four Higgs doublet models 

The most general superpotential with four higgs doublets is given by 

W = g(hiiJi + hsi/a)^' + Q(h2i/2 + ^H^W + L(h^i/i + h^i/g)^' 

+ /^12-f^l-f^2 + /^32-f^3-f^2 + l^uHiH4 + ll^4^H^H4^, (1) 

where Q stands for quark doublets, D'^ for down quark singlets, Lf^ for up quark singlets, 
L for lepton doublets, E'^ for charged lepton singlets, and hj are the Yukawa matrices 
(family indices are omitted). The Higgs doublets i^i, and H4 have hypercharges 
— 1 and +1, respectively. 

Including soft SUSY breaking terms the effective potential for the Higgs fields is 

V = rnlHlHi + rnlHlH2 + rnlHlH3 + rnlHlH4 + 
+ {ml^HiH2 + h.c.) + {m%H^H2 + h.c.) + 
+ {ml^HiH^ + h.c.) + {ml^H^H^ + h.c.) + 

+ {ml^HlH^ + h.c.) + {ml^HlH^ + h.c.) + T/^^^ + AF, (2) 
where V^^^ contains the D-terms and the radiative corrections. For the neutral 
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components (f)i of the doublets one has 

V^'''' = lig' + g")[4<Pi + 4h-4<p2-<Pi^,?. (3) 

The radiative contributions AV are generated by SUSY breaking effects. In our argu- 
ments it will suffice to consider the terms derived from large top (and possibly bottom) 



quark Yukawa interactions |T6| 



AV ^ ^J-^l,nlM^)-l]-,ntM^)-ln, (4) 

where q = t,b, m'f = |/i2t02 + Kt^'Al'^, ^1 = \hib<l)i + hb^'sl'^, and w?- = ml + m^. 

Our first comment about the viability of 4HD models should make reference to the 
size of FCNCs via Yukawa interactions. If the Yukawa matrices in (|lD are uncorrelated, 
there is no reason to expect that the unitary transformations defining mass eigenstates 
also diagonalize (in flavor space) the couplings to the extra Higgs doublets. This would 
introduce unsuppressed FCNCs at tree level. The observed pattern of quark masses and 
mixings, however, strongly suggests the possibility of an approximate flavor symmetry as 
the origin of the hierarchies required in the Yukawa matrices. In the simplest scenarios 



T7| , |18] the effect of such a symmetry would be to generate fermion matrices with off- 
diagonal elements of order 0(y^m^my/f), where rrti is the mass of the ith quark and 
V is the weak scale. If the extra Higgs doublets are a replica (with respect to this 
flavor symmetry) of the first doublet, they will introduce Yukawa matrices with the same 
approximate structure. In that case the smallness of Yukawa couplings is enough to keep 
all FCNC within the experimental limits. In particular, for extra Higgs masses around 
1 TeV the tree-level contributions to K — K and B — B mixings would be of the same 



order as the standard contributions |T8|, |T9l . In our 4HD model we will assume this type 
of approximate flavor symmetry at work. 

Our main motivation to study 4HD models concerns the origin of CP violation. In 
these models explicit CP violation seems even more inconvenient than in the MSSM, 
due to new processes mediated by the Yukawa interactions described above (the approx- 
imate flavor symmetry would not explain, for example, the small value of ex [0)- '^^^ 
possibility of SCPV in 4HD models has been addressed in a recent paper [^. There we 



assume that all the parameters in the Lagrangian are real and the CP- violating phases 
appear via VEVs ViC^^^ of the Higgs fields. We showed that at tree level {i.e., AV = 0) 
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the minimum equations for the phases can be solved in terms of a simple geometrical 
object but the remaining conditions are then uncompatible. Namely, after a redefinition 
of masses and fields that cancels the terms mf^ and m^^, the four tree-level minimum 



conditions for the VEVs Vi with nonzero phases read 

""'dv, - ^^^^^ miX4 Mv)+^^'^^^"" 

^2— = vjlml-ml^ml^hi-v) - g{v)]=Q 



' dv- 
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t;4— = + mi4m34/i(v) - 5((v)] = 0, (5) 

where ^'(v) = + g''^) [vl + ^1 — vl — and 



h(y) = \lml2mli - ml^ml2 



\ 



1 „ ..2 



(6) 



Since the four equations above depend on only two combinations of VEVs, g{\) and /i(v), 
there will be no solution [i.e., phases different from or tt) unless a fine tuned value of 
the mass parameters is imposed. Moreover, if this fine tuning were used it would imply 
the presence of two massless scalar fields. 

The effect of the radiative corrections is twofold: they relax the ammount of fine 
tuning required in the equations above, and they generate masses for the two massless 
modes. For these two effects to be sizeable we need (see AV in Eq. ^) large squark masses 
and at least two large Yukawa couplings. These could be the two top quark couplings 

and h^t or one top (/i2t) plus one bottom (hib) coupling. However, for < 5 TeV 
and Yukawas smaller than ^1.2 (as required to avoid Landau poles before the Plank 
scale) we find that the two light scalar fields have masses smaller than ^ 30 GeV. In 
consequence we conclude that in 4HD models with all the parameters real the presence 
of nontrivial complex phases in the Higgs VEVs implies two scalar fields apparently too 
light. (A more detailed examination of the parameter space might show, however, that 
this possibility is not entirely excluded by current limits on the masses of the scalar 
fields.) The situation here is then similar to the MSSM (where the allowed mass of the 
light scalar field is already excluded |lT[) or the singlet model (which relay on radiative 
effects to give mass to a mode with negative tree-level mass |T2|). 
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In 4HD scenarios there is, however, still another possibility which seems consistent 
with the idea of SCPV. It requires that the four Higgs mass parameters /i in the super- 
potential (or, equivalently, the six parameters mfj in Eq. |^) are allowed to be complex. 
This could be justified since the Higgses are the only superfields which are not protected 
of mass contributions by the gauge symmetry. They could acquire their masses in an 
intermediate scale, via complex VEVs of singlet fields with no sizeable effect on the 
rest of the low-energy effective Lagrangian. We will not assume complex phases on all 
soft SUSY-breaking terms, since in principle these singlets do not couple to gauginos or 
squarks (we will neglect the possibility of further phases or new contributions to SUSY- 
breaking parameters due to the presence of nonsinglet heavy fields 0]). The hypothesis 
of complex Higgs masses, consistent with a soft origin of CP violation, is not possible 
in the MSSM or the singlet model, since there (unlike here) all the Higgs masses can be 
made real by field redefinitions. 

In the next sections we study the implications of a 4HD model where all the param- 
eters in the initial Lagrangian are real except for the Higgs mass parameters. 

3 Definition of the model 

The approximate flavor symmetry described in the previous section suppresses all FCNC 
amplitudes to acceptable limits. However, multi-Higgs models face a potential problem 
also with CP violation: if the Yukawa couplings are complex with phases of order one, 
CP- violating signals in K physics would be too large. In particular, ex would be typ- 
ically two or three orders of magnitude larger than observed [|l9l. Thus it seems that 
a general model of many Higgs doublets requires another ingredient in addition to the 
flavor symmetry. Its effect should be either a suppression of the Yukawa couplings of 
the new doublets, or to make the complex phases small. The first approach is typical 
in models with natural flavor conservation (NFC) p2[, whereas a natural suppression of 
the phases has been obtained in the superweak model with SCPV proposed in 0. In our 
scenario these two effects will be achieved by the action of a global symmetry. 

We will assume that the effective Lagrangian of the model obeys an approximate 
Peccei-Quinn like symmetry with the following assignment of charges PB| : 

Q{H,) = +1 Q{H,) = -1 Q{D') = +1. (7) 
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All other superfields have zero charge^. The symmetry is approximate in the sense 
that couplings of operators violating the symmetry are suppressed by powers of a small 
parameter epQ. In this section we discuss the impact of this global symmetry first on the 
Higgs scalar part of the potential and then on the Yukawa sector. 

The assignment of charges tells us that in the scalar potential m^4, m|2, m^g, m\j^ are 
suppressed by epQ. ^"^^ "^34 remain unsupprcssed (of the order of the SUSY- breaking 
scale nis < 1 TcV). For easy reading we will write the suppression factors explicitly; for 
example becomes epQm\2-, where 771^2 = 0{m^). The tree-level scalar potential (we 
neglect radiative corrections in the following) involving only neutral Higgs fields is then 
given by 

V 








(8) 

As explained in the previous section, we assume that the mass parameters mf^ are 
complex. The first two mass matrices above can be diagonalized through two unitary 
transformations of order epg of the scalar fields: 

(9) 

The quartic term in the potential will not change its form and can be obtained just by 
replacing unprimed by primed fields. The relative size of the four complex masses in 
the third mass matrix above will stay the same {i.e., the off-diagonal elements are still 
suppressed by epg). A phase transformation of the fields 0j can be used to remove three 
of the four phases, leaving only one phase a in the scalar potential. Dropping the prime 
to specify transformed quantities, the mass terms read 




^ There should also be charge assignments in the lepton sector, e.g. Q{E'^) = We will comment 
on this possibility when discussing the electron EDM. 
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where the masses are now real and the phase a of 777,^4 has been written explicitly. We 
assume a to be of order one, since there is no symmetry reason for it to be suppressed. 
It is easy to see how the Higgs field redefinitions above change the mass parameters /ijj 
in the superpotential W (these parameters arc relevant since they will appear in scalar 
trilincars). We obtain ^ly^ and /i34 real up to order epg whereas /ii4 and /i23 will be mass 
coefficients with arbitrary complex phases but supprccsscd by a power of epg. 

We now go to the minimization of the Higgs potential. In particular, we want to find 
what are the relative size and the phases of the scalar VEVs suggested by the approximate 
symmetry. We write 

1 1 A 

< 01 >= ; < 03 >= -^^"^ ' (11) 

and 

where a global hypercharge transformation has been used to rotate away the phase of 
< 01 >. A detailed discussion of the minimum equations can be found in the Appendix. 
The results are the following. For nonzero values of the phase a the minimum is allways 
complex. The suppression in terms of epQ of the mass parameters determines the order 
of magnitude of the VEVs and phases: 

V\,V2 = 0{v) 
fa, t'4 = 0{epQv) 

5s, 54 = 0(1) (13) 

where v denotes the weak scale. A remarkable feature of the model is that one can 
understand its structure in terms of an expansion in epg from the model with just two 
doublets. In the limit epQ — the sectors {Hi, H2) and {H^, H4) decouple; the minimum 
gives then equations for vi and V2 identical to the VEVs in the MSSM, whereas vs — V4 — 
0. The phase 82 is then zero, while the phases S3 and 84 are irrelevant. Turning on a small 
value of epQ gives (proportionally) VEVs to the extra pair of scalars. Simultaneously it 
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allows a nonzero value of a, which translates into unsuppressed complex phases in the 
{H^^H^) sector and a phase of order epg in the {Hi,H2) standard sector. The mixings 
between the two sectors are small: either in a basis of scalar mass eigenstates or in a 
basis where f 3 = = (useful when discussing FCNCs via Yukawas), both are obtained 
from the original basis just by unitary transformations of order epQ. 

We now turn to the Yukawa sector of the theory. The charge assignments dictates 
that the matrix h2 is unsuppressed, hi and h4 are suppressed by a factor of epg, while 
is suppressed by epg. Making this suppression explicit the Yukawa sector for the quark 
fields reads 

= Q{epQhiHi + elQh3Hs)D'' + Q{h2H2 + epQh^H,)U\ (14) 

where now hj (z = 1,...,4) just carry the suppression from the flavor symmetry. In 
our model the Yukawa couplings of the initial Lagrangian are real. However, the field 
redefinitions performed to leave only one phase a in the Higgs potential will also redefine 
the Yukawa couplings and introduce complex phases. As we show, these phases translate 
into a CKM phase and complex FCNC couplings which are naturally suppressed by the 
approximate global symmetry. 

We first performed the unitary transformations in Eq. (P), which redefine the fields 
Hi by (complex) factors of order epg. They imply a redefinition of the Yukawa matrices 
which introduces phases of order epg in hi and h2 and of order one in and h^. Then 
we performed the (order one) phase redefinitions of the Higgs doublets that make all 
mass parameters real except for ml^. This translates into overall phases of order one 
multiplying the Yukawa matrices hj. However, we can still redefine the quark fields 
and absorb the phases which multiply hi and h2 (the leading Yukawa couplings). The 
net result is that the Lagrangian in ([l^) expressed in terms of the Higgs fields used to 
minimize the potential has real (up to order Cpg) couplings in hi and h2 and arbitrary 
complex phases in hs and h4. 

After spontaneous symmetry breaking, the structure of VEVs in ([131) s uggests (note 
that vi amd V2 are not suppressed by powers of epn) tan 13 = \ h^rx^ = 0(1) and 

V '"l~^'"3 

^PQ = O^mb/rrit). Thus, the approximate global symmetry is used to accommodate the 
small ratio rrih/mt, while the hierarchy between generations of the same charge is left 
to the flavor symmetry (the flavor symmetry would be exact in the limit with only the 
third generation being massive). The complex Yukawas and h4 and the corresponding 
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VEVs (f3 and ^4) are both suppressed by a power of epg, whereas the leading Yukawas 
and VEVs are real up to order epg. We then obtain quark mass matrices where all the 
entries have complex phases of order e^Q. In consequence, the complex phase in the 
CKM matrix is also of order epg. 

It is also easy to see what is the pattern of FCNC and CP violation via scalar 
exchange predicted by the model. It will be convenient to define a basis where only two 
of the four Higgs fields develop VEV and then only the second pair of scalars mediates 
FCNC processes. Again, this involves a unitary transformation of order epQ, which 
leaves almost decoupled the extra pair of Higgses (essentially {H^, H4)). In addition, the 
mixings in the scalar mass matrices between the two sectors are also suppressed. The 
overall suppression by a power of epg, when added to the one with origin in the flavor 
symmetry, renders these tree-level FCNCs smaller than CKM (box) diagrams typically 
by a factor of epg. In particular, the K — K and B — B mixings are dominated here by 
the standard contributions, like the VT-exchange box diagram in Fig. |1]. This fact will 
distinguish our scenario from typical multi-Higgs models with soft CP violation where 
the tree-level superweak interactions are the main source of flavor changing processes 0. 
On the other hand, the sector (7/3, H^^ involves arbitrary phases in Yukawa couplings 
and scalar VEVs. Although suppressed by a factor of Cpg, these couplings can be the 
dominant source oi CP violation through diagrames like the one shown in Fig. ^. In 
particular, as shown in the next section, they are the main source of complex phases in 
K physics and compete with box contributions in B physics. 

We need as well that FCNC contributions via SUSY particles (wino and gluino box 
diagrams) are within the experimental limits, which in general requires certain degree 
of squark- quark alignment and squark degeneracy. In fact, the squark- quark alignment 
could appear here as a natural consequence of the flavor symmetry ||2^] . Since we assumed 
no complex phases in soft-SUSY parameters others than Higgs masses, their CP violat- 
ing effects will follow the same pattern described above. We explore these and other 
phenomenological implications of the model in the next sections. 
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4 CP violation in K and B physics and the neutron 
EDM 



K physics. As our first example we look at the K system. In this scenario the FCNC 
processes via Yukawa interactions (see previous section) are highly suppressed, and the 
dominant contribution to Re AM12 comes from the box diagram in Fig. |I|. The leading 
imaginary contribution to AM12, however, will come from the neutral Higgs exchange 
in Fig. ^. The fiavor- changing Yukawa couplings are complex (with phases of order 1) 
and generically suppressed by the Peccei-Quinn and the fiavor symmetries (for example, 
in Fig. H the couplings are of order epQ^Jm^nllv). When the mass of the exchanged 
scalar is around 1 TeV, this (complex) diagram is roughly suppressed by Cpg with respect 
to the box diagram in Fig. [l|. Since the CKM matrix and then the box contributions 
are approximately real, the leading contribution to Im (AM12) comes from the Higgs 
exchange in Fig. 0. On dimensional grounds, the CP violating parameter (see |^| 
for definitions and notation) in the K system and epg are related: 

1 Im AM12 2 
''''''^2v^ReAMi2'^'^«- ^ ^ 

The parameter epQ sets the overall strength of Yukawa couplings of the Higgs doublets 
and suggests the order of magnitude of all the scalar VEVs. In particular (see section 
3), one expects tan/5 = 0(1) and epq = 0{mb/mt). Then the relation above establishes 
ex ~ 10~^, as experimentally required. 

Other sizeable contributions to ex may come from SUSY box diagrams with chargino 
or gluino exchange. Both of them require large SUSY contributions to AM12 (of the 
same order as the standard box diagram). Chargino box diagrams would then give 
contributions of order ex ~ 10~^epQ, whereas gluino boxes could be as large as 
ex ~ epQ {i.e., of the same order as the dominant tree-level scalar exchange). The 
factors above derive from the suppresion in Yukawa couplings or extra scalar VEVs. 
Large gluino box contributions, however, also require large left-right squark mixing 6lp = 
Tn\^/m1 ~ 10"'^ (a naive estimate would give 5lr ~ A^rn^rns ^ io~^). In addition, the 
tree- level contributions to ex can be easily enhanced P3| assuming Higgs masses ligheter 
than 1 TeV, so the clear tendency in our model is that this type of nonstandard Higgs 
exchange provides the dominant contribution to ex- 
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In contrast, the expected value for e'^ differs in principle from the standard model 
prediction. An estimate of e'j^ can be obtained from the phase to = ImAo/ReAo, where 
Ai is the decay amplitude of a K° into two pions of isospin i (see |^ for notation). In 
particular one has the experimental constraint 

to^v^l4^||6'^|<10-^ (16) 
The dominant contribution to ReAg arises from the standard penguin diagram: 

^ sin In ^Ocn + H.c. , (17) 

where 

Ocn = {sLl.T^dL) (qRlf^T-qR) (18) 

and T" are the 3-dimensional generators of SU{3). In our scenario, however, the imagi- 
nary part of this penguin diagram is suppressed by the smallness of the phase (of order 
~ 10"'^) in the CKM matrix. Since the standard model prediction for to is of order 
S13S23/S12 ~ 10^^, we obtain a first contribution of order 10"^. Other contributions to to 
may come from penguin diagrams with chargino and stop (Figure ^ and tree-level dia- 
grams with charged scalars. The first contributions have been studied in in the context 
of SUSY models with SCPV. It is found that they are typically of order to ~ 10~^sin5, 
being 6 the complex phase in the VEV of H2 (the scalar field giving mass to the top 
quark). In our model the Higgs fields H and H in Fig. |^ can correspond to H2 or H^. 
In the first case the phase 6 is of order epg, and in the second case the same degree 
of suppression comes from the small VEV and the small Yukawa coupling of H4. The 
contributions via exchange of charged Higgses have been analyzed in in the context 
of two-Higgs doublets models. In our model either the Yukawas are almost real (for the 
doublets with dominant couplings, as it happens in or the Yukawas themselves are 
suppressed. Hence, from the three types of diagrams we obtain contributions (taking 
e|>Q ^ 10~3) of order to ^ 10"^ or ^ ^ 10"^ (for IA2M0I ^ 1/22). 

Potentially larger contributions to e'^ are expected from gluino-mediated penguin di- 
agrams (Figure H). Although gluino masses in our model are real, there will be complex 
phases of order epQ (see discussion of chargino penguin above) in left-right squark mix- 
ings. This type of contributions to e'^ have been studied in detail in 0. There it is found 
that for complex phases in gluino masses of order 10~^ they could result in values as large 
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as ^ ~ 10^'^. An analogous result has been recently obtained in with contributions 
from squark mixings of CKM type and small CP- violating phases (which are natural 
in our scenario). These values are only obtained, however, when gluino box diagrams 
saturate the value of e^- Since we have assumed that the FCNCs are here dominated 
by standard box diagrams, we expect a value for e'j^ typically smaller. Modulo hadronic 
matrix element uncertainties, we estimate 

^f«10-^-10-^ (19) 
with the possibility to consistently increase this value via gluino penguin contributions. 

B physics. We consider now CP violation in B physics Although today the only 
observed CP violation is in the K system, the standard model predicts clear signals in B 
decays that should be observed in the near future. These CP asymmetries are generally 
parametrized in terms of the complex phases Ajg, which in turn depend on the product 
of phases in three amplitudes: the direct h decay, the B — B mixing, and (possibly) the 
K — K mixing. In CKM scenarios the phases Ajg are constrained by unitarity and have a 
simple geometrical interpretation (these predictions are not expected to change much in 
minimal SUSY models). However, in our 4HD scenario the three amplitudes above have 
complex phases of order epg (see below), and the predictions change to the extent that 
no CP asymmetries will be observed at the projected B factory at SLAC. 

To see why this is so, we will first consider the decay amplitude of a 6 quark into lighter 
flavors. The main contribution corresponds to a tree-level diagram with W exchange. 
Since it will be proportional to elements of the CKM matrix, its imaginary component 
will be suppressed by a factor of epg. The decay via charged Higgs are suppressed by the 
same factor due to the relative smallness of their Yukawas and the smallness of the mixing 
(in the scalar mass matrix) between the standard and the extra Higgs sectors. In non- 
SUSY models with NFC the second effect can give significant contributions (proportional 
to nit) in B decays and in flavor-changing processes ||^. 

The main contribution in this model to B — B mixing AM^-^ comes from the stan- 
dard box diagrams, and is proportional to CKM elements. The tree-level diagrams with 
exchange of neutral scalar give contributions which are suppressed by the flavor and the 
Peccei-Quinn symmetries, with a relative factor of Cpg with respect to the box diagrams. 
Both types of diagrams introduce imaginary components of order epg with respect to 
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the main real component, and the contribution to the complex phase Xiq from B — B 
mixing is negligible (of that order). As discussed above, the same conclusion applies to 
the contribution from K — K mixing. 

In consequence, in this scenario one expects that all CP asymmetries in B decays 
negligibly small (of order epg ~ 10^'^). This type of prediction is shared, for example, by 
non-SUSY multi-Higgs doublet models or SUSY models with real Yukawa matrices 

0, In CKM scenarios the situation is essentially different. There the smallness of CP 
violation in the K system is atributed to the smallness of the CKM elements involving 
the third family of quarks, whereas CP-violating asymmetries in the B system are large: 
the B — B mixing and the h decays are proportional to elements of the CKM matrix with 
arbitrary complex phases. The absence of CP asymmetries at the SLAC B factory would 
point to a non-CKM origin of CP violation, and many-Higgs doublet model (SUSY or 
non-SUSY) would appear as a natural candidate. 

Neutron EDM: As in usual SUSY scenarios, the prediction of our model for the neu- 
tron EDM dn is much larger than in the non-SUSY standard model. In the MSSM 
the explicit phases ip in SUSY-breaking gaugino mases and scalar trilinears introduce 
contributions which roughly require a suppression of 2 or 3 orders of magnitude: dn ~ 
^'-'"^^(txkF^ ) ^ [0]' These diagrams are also present in our scenario, but their con- 
tribution has a natural suppression of order respect to the MSSM. The origin of 
this factor is (again!) either the smallnes of the complex phase 62 (of order Cpg) in the 
two standard Higgs doublets, or the combined relative smallness of VEVs and Yukawa 
couplings (both suppressed by a factor of epg) of the two extra doublets. 

To illustrate this fact, let us consider the contribution from the the one- loop chargino- 
squark diagram (Figure H). When H corresponds to Hi, then the complex phase in the 
Yukawa coupling is suppressed. When H is H^, then the VEV and the Yukawa coupling 
are of order epg. In consequence, in this 4HD one expects 

dn ~ lO^^^elg e cm ^ 10'^^ e cm , (20) 



a value which is close to the present experimental limit \dn\ < 1.2 x lO^^^e cm [^]. 

Here we also comment on the lepton sector. We still have the freedom to assign a 
global symmetry charge to E'^ (or even L). For simplicity let us consider E'^ = +1, 
which would be consistent with m,- = 0(1711,). In this sector all FCNC processes via 
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nonstandard scalars will be completely negligible (the size of Yukawas suggested by the 
flavor symmetry would be enough to control all these processes). The pattern of CP 
violation will be analogous to the one discussed in the quark sector, with the relevant 
complex phases suppressed by a factor of Cpg. The leading contribution to the electron 
EDM comes from a diagram similar to the one shown in Fig. ^ If all the superparticles 
have comparable masses it is expected that de ~ 10^^c?„ so that in our model the 
electron EDM is not far from the present experimental limits. 

5 Other phenomenological implications 

As shown by Flores and Sher in |jl3|, the presence of a light neutral scalar field (with a 
tree-level mass smaller than Mz) is a prediction shared by all SUSY models with Higgs 
doublets only, regardless of the number of doublets. Since in the limit epq ^ the 
scalar sector of our 4HD model essentially coincides with the MSSM, we expect small 
corrections to the standard predictions. 

To see how these corrections arise we will first consider the model with a = and, 
consequently, with all the VEVs real. The approximate symmetry dictates that v ^ Vi ^ 
V2 and epqv ^ ^ v^. We can perform two rotations of order epq of the Higgs fields 
(one in the space 0i — 03 and another in 02 — 04) in such a way that fs = f4 = 0. It is 
then straightforward to find the mass 4x4 matrix M| 



for the CP-even scalar fields hi. The 2x2 matrix Mq corresponding to hi — /12 is 
identical to the one obtained in the MSSM, with an eigenvalue m\ smaller than M| 
and another ^ m^. The submatrix M| in the /13 — /i4 sector has two eigenvalues 
of order m^. The Peccei-Quinn symmetry forces all the elements in Mf = 0{epQm1) 
and, through mixing, tends to lower the lightest eigenvalue in M| by terms of order 
e^PQml. For nonzero values of a the scalar VEVs will be allways complex (see Section 
3), introducing mixing betweeen CP-odd and CP-even states. Due to the approximate 
symmetry, however, the mixings of the lightest scalar field with CP-odd states are small 
and introduce corrections of the same order. In consequence, we conclude that these 
corrections do not change significantly the tree-level bound irth < Mz (for epg = 10~^ 




(21) 
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and rris = 500 GeV this bound is lowered by less than 2 GeV). However, we expect 
radiative top quark effects to be much more important. 

We note that the spontaneous brealcing of the (approximate) global symmetry do 
not introduce light fields. The reason for this is that in the limit of exact symmetry 
(epQ = 0) the only VEVs breaking the symmetry (^3 and V4) go to zero too: there are 
no light pseudo-goldstone states because the size of the spontaneous and the explicit 
symmetry breaking terms is of the same order. 

It is also easy to see that this model accommodates the small ratio mh/mt without 
need of fine tuning to avoid too light charginos fS^] (in the MSSM, a small mass ratio 
irtb/mt based on a large value of tan/5 implies such fine-tuning problem). The chargino 
mass matrix is here 



( /ii2 epQiiiA ^ \ 

g^pQV 

gv_ gepQV 
\ V2 V2 / 



(22) 



where we used the VEVs in (|^) and denoted the gaugino mass by M. This structure 
has no light eigenvalues. 

Another possible implication of 4HD models concerns the value of Rb. Within the 
standard model, the partial width of the Z boson to bb seems to be very sensitive to top- 
quark radiative correction. For the top observed in CDF the predicted value is well below 
(a three-o" deviation) the present experimental limits |]3T[. In minimal SUSY scenarios 
the main correction results from the balance between Z vertices with Higgs-top and their 
SUSY partners, and the anomalous value of Rb can be alleviated for light charginos and/or 
light stop scalars B^. In 4HD SUSY models the situation is similar (especially in our 
scenario, due to the global approximate symmetry assumed in the Yukawa sector), with 
more freedom than in the MSSM to adjust the corrections. Note, for example, that large 
bottom Yukawa couplings do not imply necessarily a large value of tan/5 (= • ' — — -' 



Our last comment concerns the strong CP problem. In the model under consideration 
there are (tree-level) contributions to 6 of order epg ~ 10~^, a value much bigger than 
the present experimental limit {6 < 10^^). It seems possible, however, that the interme- 
diate scale used to break CP would also define a realistic axion scenario. Some of the 
ingredients (a Peccei-Quinn symmetry or singlet VEVs breaking the global symmetries) 
are already present in the model. Of course, for this scenario to work other requirements 
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(on the dimension of the operators breaking the anomalous Peccei-Quinn symmetry, on 
the ratio of the scales involved,...) are also needed. 

6 Conclusions 

The origin of CP nonconservation in SUSY models provides a good reason to explore 
nonminimal extensions. In the usual MSSM scenario CP-violating phases occur in two 
different sectors: in Yukawa couplings, where they would be responsible for CP violation 
in K and B physics, and in SUSY-breaking terms (gaugino masses and scalar trihnears), 
where they would induce too large fermion EDMs unless suppressed by two or three 
orders of magnitude. 

We have presented here an extension of the MSSM with four Higgs doublets where the 
complex phases appear spontaneously, induced by explicit phases in Higgs masses. An 
approximate Peccei-Quinn symmetry almost decouples the pair of extra liiggs fields, but 
their small couplings (also suppressed by the flavor symmetry) turn out to be responsible 
for all CP- violating phenomena. In particular, tree-level FCNC diagrams are irrelevant in 
ReAMi2 but responsible for ex- The resulting CKM matrix of the model has a negligible 
complex phase of order e^pq ^ 10~^. This suppression appears in all CP signals either 
from small phases in the dominant scalar sector or from small ratios of VEVs and Yukawa 
couplings in the extra sector. 

On dimensional grounds, the parameter epq specifying the violation of the Peccei- 
Quinn symmetry sets: 

• the ratio mb/rrit ~ epq, and the relative suppression of the Yukawa couplings of the 
extra Higgses {h^/hi !v h4/h2 ^ ^pq)', 

• the parameter ek ~ Cpg and the ratio e^/e^ < Cpg (with a preferred value (10~^ — 
10-^)6|,g); 

• the neutron EDM lO^^^ep^e cm, being 10~^^ an estimate for typical SUSY- 
breaking parameters; 

• and the CP violating asymmetries Xiq ~ epg involved in B physics. 

In consequence, a neutron EDM close to its present experimental limit, negligible CP- 
violating effects on B physics, and a small value of the e'^ parameter could be regarded 
as typical predictions of the model. In addition, we have estimated the effects of the 
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extra sector on the mass of the hghtest neutral scalar and commented on other aspects 
of the model {Rb and the strong 9 parameter). 

We think that 4HD models constitute an interesting possibihty in SUSY extensions 
which, however, seems almost absent in the literature. We have defined a scenario where 
CP violation is brought under control in a consistent way (due to the action of an approx- 
imate symmetry), in contrast to SUSY models where the complex phases are assumed 
small without explanation. Although wc have analyzed a particular model, we think that 
it contains essential ingredients which may be shared by any satisfactory multi-Higgs 
SUSY model. In a generic multi-Higgs model hard (CKM-like) CP violation seems to 
imply too large imaginary FCNCs mediated by the extra Higgs fields. This fact strongly 
suggests a soft origin of CP violation. Then the problem of containing simultaneously 
FCNC and too large CP violation forces these models to have, for example, unobservable 
CP asymmetries in B decays, a prediction that will be tested in the near future. 
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Appendix: Solutions to minimum equations 

The vacuum expectation value of the scalar potential is 

^T/^ ^22,122, -'^22,122, 2 e, 

< V > = -rriiV^ + -777,2^2 + 2"^3'^3 + 2"^^"^^ mi2VlV2 COS 02 + 

+ ^PQml^viVi cos((54 + a) + epQml2VzV2 cos{5^ + 82) + ml^v^Vi cos{5^ -\- 5^) -\- 
+ ^{9' + 9")[vf + vl-vl-vir. (23) 

The conditions at the minimum are 

dV 

vi- — = mlvf + ml2ViV2Cos62 + €pQml^viViCos{6i + a) + vlg{v) = , 
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dV 

V2-^ — = mlv2 + ml2ViV2CosS2 + epQrnl2V3V2Cos{S3 + 62) - vlg{-v) = 



2 



dV 

^^dv' ^ mlvl + epQml2V3V2Cos{S3 + 62) + ml^V3ViCos{53 + 5i) + vlg{-v) = , 
dV 

V4T— = mlvj + ml^V3Vi 008(63 + ^4) + epQ"^?4^i^4 003(^4 + a) - vlgiy) = (24) 



dv4 

where ^(v) = |(^^ + g'"^) [vf + vl - - vj], and 
dV 



352 
dV 

dV 



mi2'"i'"2 sin 62 + epgmgaf 3f 2 sin ((53 + ^2) = , 

sin((53 + ^2) + "^34^^3^^4 sin(53 + ^4) = , 
r.r rnl^V3V4 sin(53 + ^4) + epgrnl^viVi sin(54 + a) = 0. (25) 

0*04 

First, we can estimate the relative sizes of VEVs. The solutions to (^4]) are consistent 
with either one of the pairs of VEVs (fi,t'2) or (f3,f4) being suppressed by order epg 
with respect to the weak scale. Depending on the sizes of the unsupressed parameters 
m\,m\,m\^ and rn^^rn^^rn^^ the absolute minimum will prefer one of the above choices. 
This can easily be seen from the following. Imagine for a moment that epQ = 0. The 
equations in (El) reduce to two pairs of equations, with first pair depending on the 



ratio ^ and giy), and the second on ^ and the same function giy). Thus, one pair of 
VEVs is forced to be zero. Turning back on a small epg, the terms in the Lagrangian 
suppressed by epg can make the previously trivial pair nonzero, but suppressed. The 
only question left is which pair of the VEVs is small, and this will depend on the choice of 
the unsuppressed mass parameters in the potential. We will assume that the parameters 
are such that {vi,V2) are unsuppressed (and of the order weak scale), while (f3,f4) are 
of the order epg times the weak scale. This assumption does not involve fine tuning but 



only "halves" the available parameter space. For the phases 5j, it was shown in pO[ that 
in the limit a = there is no nontrivial solution [i.e., 6i different from or vr). It is easy 
to see, however, that for a ^ the equations (^) do not have this trivial solution, and 
complex phases are guaranteed. In particular, for a = 0(1) from the two first equations 
in (^) it follows that 62 is of order epg (modulo vr, depending on the sign of m^g) 
63 and ^4 are unsuppressed. 

In summary, the structure of values of VEVs and their phases is 

Vi,V2 = 0(v) 
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V3,V4: = 0{epQv) 
63,6, = 0(1) 



(26) 



where v denotes the weak scale. 

We now explore this structure in more detail (at first order in epg). The first equation 
in (^) gives ^2 to be of order tpg up to^ a factor of vr. Such 62 does not contribute to 
leading order (its contributions are O(epQ)) to the minimum of the scalar potential (|23|) , 
and can be neglected in the rest of equations. The minimum equations (|2^) to leading 
order are 



dV 
dvi 
dV 

dv2 
dV 

dvz 
dV 

where go(y) 



V2 



V3 



2 2 I 2 
2 + Tn^2'^iV2 



2„,2 



0, 
0, 



ninV 



3^3 







epQm^2^3t'2 cos ^3 + ml^v^Vi cos(53 + ^4) + f 35'o(v) 
m\vl + ml^v^VA cos(53 + ^4) + epQm\^viVA cos((54 + a) - vlgoiy) 

-epQml2V3V2 sin ^3 + ml^VsVi sm{5s + S4) 



,(27) 



dV 
dV 







"^34'^3t^4 sin((53 + 64,) + epQm\^viV4 sin((54 + a) = 0. 



(2^ 



The first two equations in ( pT]) are just the equations of the MSSM, and they fix vi and 
V2 in the usual way. Thus, we expect both vi and V2 to be of the order of weak scale and 
tan/5 = 0(1) {i.e. no supression by epg, and no fine tuning producing large tan/3). 

The goal now is to find the phases ^3 and 5 4 in terms of quantities c = epQW^32'y3^2 5 
/ = ml^v^Vi and y = —epgrn'^ViV^ and the angle a in order to eliminate them in the 



third and fourth equation of VEVs in ([271 ). We note that the three quantities c, /, and y 
are of order O(epQ), and so we expect ^3 and 64 unsupressed. In order to find ^3 and ^4 



we use a geometrical interpretation similar to the one devised in pO[. It is possible to see 
that the two equations ( pSD define one of the objects shown in Figure ^ (which one it is 

^ In the following we choose positive without loss of generality and thus ^ tt + ©(epg). 
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will depend whether 1/c, 1/f and 1/y can form a triangle or not). The quantities p and q 
there are not independent, and can be expressed in terms of c, /, ?/ and a. The difference 
between the two types of solutions can be understood in the limit a — 0, where only 
the trivial solutions ^3 = and ^4 = vr exist. For a = the object in Fig. |^(b) implies 
nonzero ^3 and ^4, a type of solution which requires fine tuning between mass parameters 
once it is substituted in the equations for t>3 and ^4. In consequence, for small values of 
a only the solutions of the type in Fig. ^(a) appear. When a is nonzero the fine tuning 
is lifted, and both types of objects define possible solutions to the minimum equations. 

We performed numerical solutions to the above equations when a 7^ and large (order 
1) and we found that the minima satisfy the structure given in (^). For simplicity and 
to illustrate the discussion above we will give the equations for and at first order 
in a. 

Fig. ^(a): When a ^ we see that ^3 ^ and ^4 vr, while > 1/c+l// and 
1/q — 1//- From the figure we first find cosines of the relevant angles {63 + ^4, 

63 and 64^ + a) to leading order in a. Then we are in the position to find ^3 and V4 by 
substituting these expressions into the two last equations of (p7|) 

VsZ = mlvl-c[l-^{^-^r]-f[l-^{^-^r]+vlgi.) = 0, 



fey f ' 



y 



= mlvl- f[l-^i I _J ]+y[l- I _, r]-vlgi.^) = 0(29) 

^ fey fey 

where, again, c = epQm|2f3'y2, / = 'fn'i/iV^Vi and y = —epgrn'^viv^. These equations, 
although still complicated, can be solved in and V4, (remember that Vi and V2 are 
already fixed). Then we can go back and find ^3 and ^4, thus completing the search for 
the first case. 

Fig. ^(b): In this case 1/p 1/y and 1/q 1/c, while ^3 and 64 tend to go to 
angles in the triangle with sides 1/c, 1/f and 1/y (we denote this (order 0(1)) asymptotic 
angles as ^3 and 62)- We can again find the relevant angles to leading order in a and 
then find V3 and V4 by substituting these expressions into the two last equations of ( P7| ) 

2r 2 , "^32"^34'^2 , 0/ NT , o x ^^(^s + ^4) 

= ,H,nl + - /(V)] + 2a/ f 'g ; 'I ^ . 

0V4 V2 tandftandl 

(30) 
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Remembering that vi and V2 are already fixed in terms of m\ and m\2i we see 
clearly the fine tuning for vanishing the terms in square brackets would be forced 

to vanish, implying two relations between mass parameters. However, once we include 
a 7^ these degeneracies get lifted. 
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Figure 1: Leading contribution to Re AM12. 
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Figure 3: Chargino contribution to e'^ 
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Figure 4: Gluino contribution to e'^ 
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Figure 5: Chargino contribution to the neutron EDM. 
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Figure 6: The two possible geometrical objects which represent the CP nontrivial 
solution of equations (29) generated by a nonzero soft phase a. Each object consists of 
two triangles, ABC and ADE. The sides of the triangles are AB = 1/p, BC = 1//, AC 
= 1/c, AD = 1/y, DE = 1//, AE = 1/q. (a) The object is such that the sides 1/c, 1/f 
and 1/y cannot form a triangle, (b) The object is such that the sides 1/c, 1/f and 1/y 
can form a triangle. 
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